Observation of Asymmetric Transport in Structures with Active Nonlinearities 
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A mechanism for asymmetric transport based on the interplay between the fundamental symme- 
tries of parity (V) and time (T) with nonlinearity is presented. We experimentally demonstrate and 
theoretically analyze the phenomenon using a pair of coupled van der Pol oscillators, as a reference 
system, one with anharmonic gain and the other with complementary anharmonic loss; connected 
to two transmission lines. An increase of the gain/loss strength or the number of 'PT-symmetric 
nonlinear dimers in a chain, can increase both the asymmetry and transmittance intensities. 
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Directed transport is at the heart of many fundamental 
problems in physics. Furthermore it is of great relevance 
to engineering where the challenge is to design on-chip 
integrated devices that control energy and/or mass flows 
in different spatial directions. Along these lines, the cre- 
ation of novel classes of integrated photonic, electronic, 
acoustic or thermal diodes is of great interest. Unidi- 
rectional elements constitute the basic building blocks 
for a variety of transport-based devices such as rectifiers, 
pumps, molecular switches and transistors. 

The idea was originally implemented in the electron- 
ics framework, with the construction of electrical diodes 
that were able to rectify the current flux. This significant 
revolution motivated researchers to investigate the possi- 
bility of implementing this idea of " diode action" to other 
areas. For example, a proposal for the creation of a ther- 
mal diode, capable of transmitting heat asymmetrically 
between two temperature sources, was suggested in Ref. 
PQ. Another domain of application was the propagation 
of acoustic pulses in granular systems [2J. 

A related issue concerns the possibility of devising an 
optical diode which transmits light differently along op- 
posite propagation directions. Currently, such unidirec- 
tional elements rely almost exclusively on the Faraday 
effect, where external magnetic fields are used to break 
space-time symmetry. Generally this requires materi- 
als with appreciable Verdet constants and/or large size 
non-reciprocal devices - typically not compatible with 
on-chip integration schemes or light-emitting wafers [3J. 
To address these problems, alternative proposals for the 
creation of optical diodes have been suggested recently. 
Examples include optical diodes based on second har- 
monic generation in asymmetric waveguides [3] and non- 
linear photonic crystals 5 , photonic quasi-crystals and 
molecules [6], or asymmetric nonlinear structures [7]. 
Most of these schemes, however, suffer from serious draw- 
backs making them unsuitable for commercial or small- 
scale applications. Relatively large physical sizes are of- 
ten needed while absorption or direct reflection dramat- 
ically affects the functionality leading to an inadequate 
balance between figures of merit and optical intensities. 
In other cases, cumbersome structural designs are neces- 
sary to provide structural asymmetry, or the transmitted 
signal has different characteristics than the incident one. 



In this Letter we demonstrate, experimentally and the- 
oretically, a mechanism for asymmetric transport ex- 
ploiting the co-existence of active elements with distinc- 
tive features of nonlinear dynamical systems, such as 
amplitude-dependent resonances. As a reference model 
we will use coupled nonlinear electronic Van der Pol 
(VDP) oscillators with anharmonic parts consisting of 
complementary amplifiers (gain) and dissipative conduc- 
tors (loss) combined to preserve parity-time (PT) sym- 
metry. 'PT-symmetric structures were inspired by quan- 
tum field theories [8] ; their technological importance was 
first recognized in the framework of optics 9 , where sev- 
eral intriguing features were found [9-20 . For example, 
the theoretical proposal of Refs. (T9J [20] suggested using 
nonlinearities to induce asymmetric transport. Very re- 
cently the idea of creating PT-symmetric devices within 
the electronics framework was proposed and experimen- 
tally demonstrated in Ref. [21] . PT-electronics pro- 
vides a platform for detailed scrutiny of many new con- 
cepts within a framework of easily accessible experimen- 
tal configurations [2"TH2"5] . Despite all this activity, the 
PT-symmetric Hamiltonians introduced in quantum field 
theory, optics, and electronics have been restricted to 
conservative anharmonic constituents (if any) with the 
matched gain and loss exclusively linear (see however the 
theoretical works [24]). 

Here, using the PT-electronics framework, we demon- 
strate with experiment, simulations, and theory, asym- 
metric transport from PT-symmetric structures that 
belong to this relatively unexplored class of nonlinear 
systems whose anharmonic parts includes the mutually 
matched gain and loss. Such systems, natural to the 
realm of electronics for which the van der Pol model was 
originally proposed 25], can also be implemented in op- 
tics using, for example, optical amplifiers, saturable ab- 
sorbers [35] and two-photon losses to realize the nonlinear 
PT-symmetry. 

An ideal VDP oscillator has a linear anti-damping 
at low amplitudes which is subsequently overtaken by 
a cubic dissipation at high amplitudes. In electron- 
ics this is an LC oscillator in parallel with a voltage 
dependent conductance characterized by the I-V curve 
I(V) — —V/R + bV 3 . A negative impedance converter 
(NIC) (upper left of Fig. [T]) generates a —l/R term, 
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FIG. 1: (Left) Gain and loss circuits of the van der Pol VT- 
symmetric dimer. The elements shown above are incorpo- 
rated as parallel conductances in the capacitively coupled LC 
resonators making up the dimer. (Right, color online) Experi- 
mental I-V response (circles) for the gain (red) and loss (blue) 
elements along with the corresponding NGSPICE simulations 
(solid), taken at a frequency of 30kHz, typical of the active 
range of the VDP dimer. 



and we approximate the cubic turn-around with paral- 
lel back-to-back diodes moderated by a resistor. The 
time-reversed conductance is constructed with the resis- 
tor R and the diode combination interchanged. The re- 
sulting "gain" and "loss" nonlinear conductances refer to 
their low amplitude character. The respective nonlinear 
I-V curves are shown in the right panel of Fig. [T] The 
0.36i? in series with the diodes optimizes the match to 
a cubic nonlinearity. It is important to note that only 
the parameter R is used to set the gain/loss parameter 
7 = Rr 1 L I C , while the diode turn-on characteristics 
are fixed. When comparisons are made to theoretical 
models, the voltage scaling will consequently depend on 
7- 

The schematic of the dimer circuit is shown in the 
lower left of Fig. [I] The coupled LC heart of the cir- 
cuit is identical to that used in a previous work |21j with 
LM356 op-amps serving in the NICs and signal buffers. 
The two VDP oscillators are capacitively coupled by C c . 
In principle, more complicated geometries such as chains 
of active nonlinear VT- symmetric dimers could be con- 
structed by capacitively coupling additional dimers into 
the desired topology. The experimental measurements of 
this work are limited to the dimer as a check on SPICE 
simulations used for larger systems. 

Transmission lines (TL) with impedance Zq are at- 
tached to the left (lossy) and the right (gain) LC nodes 
of the dimer to complete the scattering system used to 
perform our transport measurements. Experimentally, 
these take the form of resistances Rq = Zq in series with 
independent voltage sources, here HP3325A synthesizers, 
on the right and left sides. The incoming and outgoing 
traveling wave components associated with a particular 
TL are deduced from the complex voltages on both ends 
of Rq, as sampled by a Tektronix DPO2014 oscilloscope. 
For example, on the left (lossy) side, with Vlc the volt- 
age amplitude on the LC circuit node, and Vq the voltage 



amplitude on the synthesizer side of the coupling resistor 
i?o, the incident wave on the dimer has a voltage ampli- 
tude V£ = V /2 and the outgoing wave has a voltage 
amplitude V£ = Vlc ~~ Vo/2. An equivalent relation 
holds for the right TL terminal with the ± superscripts 
interchanged, since they refer to right or left wave trav- 
eling direction regardless of the terminal orientation. 

The scattering measurements are performed for fixed 
incoming wave amplitude set by Vq of the signal generator 
on either the left or the right side with the other side set 
to zero. The generator frequency is stepped (up or down), 
and the three relevant waveforms, Vq and Vlc on the left 
and right are simultaneously captured (the Vq channel 
on the transmitted side is zero). Harmonic components 
of each wave constituent can be independently analyzed 
for magnitude, and relative phase. Instrumentation noise 
and sample time determine the accuracy of this analysis, 
which was found to be < 1%. 

Circuit behavior was numerically modeled by the 
NGSPICE simulator [37]. Modern simulation goes be- 
yond a lumped element approach, in which circuits can 
only be modeled by idealized passive component approx- 
imations. Rather, behavioral modeling focuses on the 
mathematical relations between inputs, parameters, and 
outputs of a complicated nonlinear component; often by 
incorporating virtual dependent sources into the result- 
ing model. This allows for faster speed and higher accu- 
racy in calculations of circuit behavior. Most manufac- 
turers, in fact, provide ready-made behavioral models for 
their devices. Here, for the LF356 op-amp we heavily uti- 
lized Texas Instruments' modified Boyle model, as well as 
the parametrized diode model of Fairchild's 1N914 [28] , 

In Fig(T] circuit analysis was done in the time-domain 
for each individual element (gain or loss circuit) of the 
nonlinear "PT-symmetric dimer. Using initial DC oper- 
ating conditions, a frequency-dependent source drives a 
circuit through a transient regime into steady-state op- 
eration, at which point the voltages and currents are 
recorded. The measured I-V response curve (see Fig. 
[IJ is reproduced by the simulations. 

For the transmittances T(i>) versus the driving fre- 
quency v in Fig[2j the steady-state time-domain simula- 
tions of the 'P7"-symmetric dimer in the scattering con- 
figuration are similarly obtained. Fourier analysis is used 
to extract the relevant frequency-dependent voltages and 
currents, which are then used to calculate the scattering 
parameters. As a check of the accuracy of our numeri- 
cal approach, we have also extracted the transmittance 
using a nonlinear harmonic balance circuit analysis (29) . 
We have confirmed that the results are identical, within 
numerical accuracy, to the ones obtained from the time- 
domain analysis. In Fig. [2]we also report the experimen- 
tal left and right transmittances for the 'PT-symmctric 
VDP dimer. The overall shape of the measured trans- 
mittances reasonably match the numerical simulations. 
The deviations are likely associated with a small para- 
sitic inductive coupling. 

A striking feature of the results of Fig. [2]is the fact that 
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FIG. 2: (Color online) Experimental data of transmittances 
Tl and Tr for a left (solid blue line) and right (solid red line) 
incident waves. The results of the numerical simulations with 
NGSPICE for Tl(Tr) transmittances are shown as dashed 
blue (red) lines. Upper inset: The ratio between the trans- 
mittance associated with the third harmonic and the 

fundamental frequency T^ 1 '. Here, rj = 0.031 and 7 = 0.15. 
Lower inset: Experimental T(u) for the same r\ but smaller 
gain/loss parameter 7 = 0.11. Note that as 7 decreases, the 
asymmetry and the transmitted intensity are reduced. 



the transmittance from left to right Tl{v) differs from 
the transmittance from right to left T R {v), i.e. Tl 7^ Tr. 
The phenomenon is more pronounced closer to the res- 
onance frequencies of the corresponding linear problem 
and is the main result of this paper. This asymmetry is 
forbidden by the reciprocity theorem in the case of lin- 
ear, time-reversal symmetric systems. In fact, it is not 
present even in the case of linear 'PT-symmetric struc- 
tures [22 , although the time-reversal symmetry in this 
case is also broken. At the same time, a Hamiltonian 
nonlinear medium by itself cannot generate such trans- 
port asymmetries. Furthermore, we find that increasing 
the gain/loss parameter 7 which is responsible for the 
asymmetric transport, maintains or even enhances the 
transmitted intensities (compare the lower inset of Fig. 
[2] with the main panel). This has to be contrasted with 
other proposals of asymmetric transport (see for exam- 
ple Ref. [7]) where increase of asymmetry always leads 
to reduced transmittances. 

We have also developed a theoretical understanding of 



asymmetric transport by restricting our analysis to the 
basic harmonic. This approximation is justified by our 
measurements which indicate that the main part of the 
transmitted power is concentrated in the fundamental 
harmonic. An example is shown in the upper inset of Fig. 
[2]where we report the ratio T^ /T^ between transmit- 
tances of the third harmonic to the fundamental. Even 
harmonics are absent in the transmission spectra due to 
the nature of VDP anharmonicity, while for higher har- 
monics J 1 ( n>3 ) the experimental values of TW/TW are 
below the noise level of our measurements. 

Application of the first and second Kirchoff's laws 
at the TL-dimer contacts allow us to find the cur- 
rent/voltage wave amplitudes I,V at the left (L) and 
right (R) contact. We get 



dX L 

(TLr 
dr 



7(1 



vi)^ + v i + (i + c) 



d 2 V L 

dT 2 

v ft x \ d2VR 



d 2 V R 
C ^' (1) 
^ d 2 V L 



where the dimensionless current /voltage amplitudes X, V 
at the lead-dimer contacts are defined as Il/r = 



V = VnV 



QVL/R- 



The dimensionless time is r 




and r\ = \/L/C /Z is the dimensionless TL conductance, 
while we have also introduced the dimensionless capaci- 
tance c = C c /C, a measure of the intra-dimer coupling. 

At any point along a TL, the current and voltage deter- 
mine the amplitudes of the right and left traveling wave 
components. The forward V^/r anc ^ backward V~[/r wave 
amplitudes, and the voltage Vl/r and current Il/r at 
the TL-dimer contacts satisfy the continuity relation 



Vr,o)e- Mr +cc; 



Vl/R - [ V L/R T- "L/R 

1l/r = (vt /R -VZ /R )e-*" + cc 



(2) 



Note that Eqs. ([I]) contains nonlinear terms on the 
r.h.s. which are responsible for harmonic generation. 
However, as suggested by the experimental data, we can 
neglect these higher harmonics and thus restrict our an- 
alytic study to the fundamental. Keeping this in mind 
when substituting Eqs. ^ into Eqs. ([I]), we get 



(-l)^(-V+ /fl + V^ /R ) 



w a (i + c) + *(-ir + Wi - K,r + vi /R \ z )\ {vl /R + v~ /R ) + ^v+j- 



t^V(K/L ) = (-1)^7(1 - \K/L H - (1 + + 1 K/L + ^Vl/R + V I/r) 



;+/-|2 



(3) 



where we used the compact notation V R j L for V R and 

V£~ . The exponent s above takes the values s — 0, 1 for 
L,R current amplitudes respectively. 

The way to solve Eqs. ([3| is to introduce the back- 



ward transfer map [30] . The latter uses the output am- 
plitude V ou t = V R (V^) as an initial condition together 
with the boundary conditions = (V^ = 0) for a left 
(right) incoming wave. Iterating backwards, we calculate 
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FIG. 3: (Color online) (a) Transmission curves for a nonlin- 
ear PT-symmetric VDP dimer with 7 = 0.151, n = 0.0305, 
u) — 0.81, and c = 0.267. The linear PT-symmetric case 
(black line) is reported also for comparison together with the 
nonlinear all-gain (red solid line) and all-loss (blue solid line) 
VDP dimer where reciprocal transport Tl = Tr is found; (b) 
Corresponding transmittances for the cases of subfigure (a). 



1 



FIG. 4: (Color online) Density plots of the Q-factor Eq. S 
for increasing 7 (from left to right) and n (top to bottom). 
The n, 7 values are indicated in the subfigures while c = 0.267. 



the corresponding incident Vi n = (V^) and reflected 
V re fl = (Vr) amplitudes for a left (right) incident 
wave. Representative (|Vi n | 2 , |V ou t| 2 ) curves are shown 
in Fig. [3^. The associated transmittances are defined as 
T = |V ou t/Mn| 2 • Straightforward algebra gives: 

2wqc 2 

Tl= m + "(7(1 - l^l 2 ) + I ^(1 + c)) + (cu)* 

(4) 

where a = (77 - 7 (1 - |V OMt | 2 ) + ± - iw(l + c)). The 
transmittance Tr for a right incident wave is given by 
the same expression as Eq. Q with the substitution of 
7 — > —7 i.e. Ik (7) = Tl(— 7) ^ Tl(j). The origin of the 
latter inequality is due to the fact that nonlinear reso- 
nances are detuned differently for left and right incident 
waves, as seen from Eq. Q (compare black dashed and 
solid lines in Fig. ^jp). The case 7 = corresponds to a 
linear passive dimer for which = Tr in agreement with 
Eq. Q. To further highlight the importance of the in- 
terplay between nonlinearity and PT-symmetry, we also 
report in Fig. [3] the transmittances for a nonlinear VDP 
dimer with both elements having gain or loss and for a 
linear PT-symmetric dimer. In contrast to the nonlinear 
PT-symmetric structure, the transmission is symmetric. 

Since the phenomenon is nonlinear the asymmetry de- 
pends on both frequency and amplitude. To quantify its 
efficiency, in Fig. [4] we report the rectification factor 

^ _ T L (V out ) — T R (V out ) 

w ~ r L (Vout) + r fl (Vout) ( ) 

which approaches ±1 for maximal asymmetry. Note that 
increasing 7,77 broadens the regions in which |Q| is rel- 
atively large. Some representative experimental rectifi- 
cation factors Q for a single dimer N = 1 and two dif- 
ferent values of 7 are shown in Fig. [5] together with the 
SPICE simulations. The measurements and the simula- 
tions compare nicely with one another and are in agree- 
ment with the theoretical calculations shown in Fig. [4] 




FIG. 5: (Color online) (Upper) Experimental (solid lines) and 
SPICE (symbols) results for the rectification factor Q Eq. |5| 
for one (N = 1) and two (N = 2) coupled PT-symmetric 
VDP dimers for two different 7 values and 77 = 0.03. (Lower) 
The SPICE simulations for Q for moderate chains of VDP 
dimers with N = 1,2,5 (we note that for v ~ 38.5kHz the 
simulations were not converging for N = 5). As N or 7 
increase the Q-factor increases as well. 

Furthermore, increasing the number N of nonlinear PT- 
symmetric dimers leads to a considerable increase of the 
Q-factor. This can be seen in Fig. [5] even for the case of 
N = 2. Although this behavior can be formally analyzed 
by evaluating the N th order backward transfer map, a 
qualitative explanation can be provided based on our un- 
derstanding of the N = 1 dimer. Specifically, one expects 
that the increase of the number of PT-symmetric dimers 
introduces 2N nonlinear resonances (which are detuned 
differently for left and right incoming waves) , thus broad- 
ening the frequency domain for which asymmetric trans- 
port is present. 

Conclusions - Using coupled PT-symmetric VDP os- 
cillators, we have demonstrated experimentally and the- 
oretically a mechanism based on the co-existence of non- 
linearity and PT-symmetry which leads to asymmetric 
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wave transport. Unlike the linear "PT or passive nonlin- 
ear cases, this is achieved concurrent with a significant 
level of transmittance. These findings can be applied to 
photonic or phononic systems as a guide for optimizing 
design in order to achieve maximal asymmetry. 
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